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Abstract. With the large amounts of spectroscopic data available today and the very large surveys to 
come (e.g. Gaia), the need for automatic data analysis software is unquestionable. We thus developed an 
automatic spectra analysis program for the determination of stellar parameters: radial velocity, effective 
temperature, surface gravity, micro-turbulence, metallicity and the elemental abundances of the elements 
present in the spectral range. Target stars for this software should include all types of stars. The analysis 
method relies on a line by line comparison of the spectrum of a target star to a lib rary of sy nthetic spectra. 
The id ea is built o n the experi ence acquire d in developing the TGME T (|Katz^ et al.i (| 19981 ) and lSoubiran et al] 
(|2003l ')'l ETOILE l|Kat jbOOll ) and Abbo l| Bonifacio fc CaffaullioO^ I software.The method is presented and 
the performances are illustrated with GIRAFFE-like simulated spectra with high resolution (R = 25000), 
with high and low signal to noise ratios (down to SNR= 30). These spectra should be close to what could 
be targeted by the Gaia-ESO Survey (GCDS). 
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1 Introduction 

One of the major applications of spectroscopy is the determination of stellar parameters like the radial velocity 
(Vr), effective temperature (TeflF), surface gravity (log(g)) , micro-turbulence (^) , metallicity ([Fe/H]) and 
chemical abundances ([X/H]). The present and future large spectroscopic surveys are going to significantly 
increase the number of spectroscopic data to be analysed. A few examples are the Gaia-ESO Survey with 
about 160'000 st ars to be observed, Gaia with about 2'000'000 stars to be analysed for chemical abundance s 
( Katz et al.ll2004 >). RAVE with some 400'000 stars observed so far (jBoeche et all (|201l[ ) and lSiebert et all (|201l[ )) 



etc. The space mission Gaia will provide the largest survey ever, and that, in the decade to come. To analyse 
these quantities of data, automatic spect ra analysis s oftwa re is needed. Different fam ilies of softw are exist. A 
few examples a re software lik e TGMET (}Katz et al.' ('1998') and 'SoubiraneLalJ (|2003l ). ETOILE (|Kat j|200ll) . 



MATISSE (Rcc io-Blanco et al . 2006), Abbo ({Bonifacio & Caffau 2003) etc... The work presented is about the 
development of a new automatic stellar spectra anaylysis software. In its first version the software will be 
optimised for medium-resolution Giraffe spectra (VLT) and will thus be tested on Giraffe like spectra. The 
software is called SPADES (Stellar PArameters DEtermination Software) and is coded with Java. 

2 SPADES 

2.1 General idea 

The software is based on the comparison between observed spectra and a grid of synthetic spectra (with known 
parameters). Contrary to many existing software, in SPADES, the comparison between spectra is not made 
over all the spectrum but around pre-selected lines. Another particularity is that the determination of the 
stellar parameters does not use any equivalent width mesures but is based on profile fitting methods. Another 
important characteristic of SPADES is that it determines elemental abundances. The general idea is as follows: 

For each parameter to be determined, one or several methods of determination (diagnostics) are available. 
One is chosen by the user. The list of diagnostics for each parameter are: 
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• Vr: cross-correlation in direct space with a template 

• TefF: excitation equilibrium or Balmer lines profile fitting. 

• log(g): ionization equilibrium or strong lines (e.g. Mgib green triplet) profile fitting. 

• [Fe/H]: Fe lines profile fitting. 

• [X/H] : X lines profile fitting. 

• ^: empiric calibration or nulling the AW — /(reduced equivalent widths) function slope. AW being the 
residuals of the difference between the observed and synthetic line. 

The diagnostics so far implemented and tested are detailed in the next section. Each of the parameters, 
Teff, log(g), [Fe/H] and ^ atmospheric parameters is determined one by one assuming all the others known. 
The process is iterated until convergence. The elemental chemical abundances are then determined. 



2.2 Diagnostics 

2.2.1 Effective temperature : TefF 
Ha wings fit method 

The first step is defining the reference grid to be used: a ID (in the parameters space) reference spectra grid 
is defined. This grid varies over Teff only, the other parameters beeing fixed to their input values. This grid 
is read or calculated by interpolation based on a pre-calculated reference spectra grid. The analysis is limited 
to the Ha line, more generarely to a spectral range around this line. This range will be now called the "Ha 
spectral domain". For each reference spectrum (TefF value), the Ha spectral domain continuum is fitted to the 
studied Ha spectral domain continuum. An example of superposed continuum fitted Ha spectral domains is 
presented in Fig. [T] Left. The spectral domains differ by their TefF. 

Halpha 

continuum adjusted and normalised to 1 spectral domain Tefl by Halpha wing fit function 

Observed versus reference 
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Wavelengtlis{A) Teff(K) 



Fig. 1. Left: The reference spectrum Ha spectral domain in red is superposed to the studied spectrum Ha spectral 
domain in black. The continuum pixels are in blue in Fig. [T]. The Ha wings pixels are in green. Right: The residuals 
as function of the reference spectra Teff. The result Teff is the value that nulls this function. 

The wavelengths ranges used for the continuum fit are pre-defined. The corresponding pixels are in blue. 
The green pixels represent the Ha wings. The wavelength limits of the wings are also pre-defined. Over the Ha 
wings pixels the residual is calculated as such: 

That is done for all the spectra in the ID reference spectra grid, s = f{Teffref) is thus constructed. An 
example of this function is presented in Fig.[T] Right. 
The result Teff is the one that nulls this function. 
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Excitation equilibrium 

In this method a hst of pre-sclccted Fel hnes is used. As in the previous method, a ID (in the parameters space) 
reference spectra grid is defined with Teff values varying around the input value while the other parameters 
are fixed to their input values. Around each central wavelength a spectral domain is cut in the reference and 
observed spectra. The line and continuum limits are determined automatically. For each reference spectrum, 
the spectral domains around the pre-selected lines are continuum fitted to the continuum of their corresponding 
spectral domain in the studied spectrum. An example of superposed, continuum fitted spectral domains is given 
in Fig. m 



Continuum normalised to 1 and adjusted spectral domains: 
Observed versus reference 
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Fig. 2. In red the spectral domain of a Fel line in the reference spectrum superposed to the studied spectral domain of 
the same line, in black. The studied and the reference spectra are at different Teffs which explains the difference between 
the line. The continuum pixels are in blue. The line pixels are in green. 

For each reference spectral domain, and over the line pixels (green pixels) the residual is measured as such: 

= Epixels -i^obs - iref) 

This measurement is done for all the used lines. Let AWn be the residual of the nth line. For each reference 
spectrum, the lines AWn arc plotted as a function of the lines respective excitation potentials An example 
of this function is given in Fig. |3]Left 

For each reference spectrum (Teff value) the slope a of this function is measured {AWn = a *?« + &)• Another 
function is then constructed: the a = f{Teff) function. An example of this function is presented in Fig. [3] 
Right. 

The result Teff is the one that nulls the a ~ f{Teff) function. 

2.2.2 Gravity: log(g) 

lonisation equilibrium 

In this method the measurement made over the Fel and the Fell lines to be used are: 

AWpel = T,pixels -(^obs " iref) and AWpell = Eptxels ^i^obs ~ iref) 
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Residuals as function of excitation potential Teff by excitation equilibrium function 
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Excitation potential Teff(K) 



Fig. 3. Left: Line residuals (AWn) as function of excitation potentials ^„ for the reference spectrum with Teff — 5800 
K. Right : The slopes a such as AWn = a * + & as a function of the reference spectra Teff. The result Teff is the 
value that nulls this function. 



The diagnostic to be analysed is : 

A - AWpel - AWpell 



where AWpei (respectively AWpeii) are the mean of the AWpei^s (respectively the AWFeii„s)ior all the 
Fel (respectively the Fell) lines. The result log(g) nulls the A as a function of the reference spectra logg 
function. 



2.2.3 Metallicity and elemental abundances: [X/H] 
Profile fit : 

In this method the measurement over the n lines of the X element to be used is : = "l^pixeis g 

The diagnostic to be analysed over all the X element lines is: = X)/mes 'l2pixeis ^'^"'""■T'''^^ 
The result [X/H] is the one that nulls the as function of the reference spectra [X/H] function. 



2.3 Tests and perfomances 

SPADES was tested by Monte-Carlo over synthetic spectra with resolution R — 25000, effective temperature 
Teff = SSOOiiT, gravity log{g) = 4.40, metallicity [Fe/H] = —1.0, with individual abundances of Ca and Ni 
: [Ca/H] = 0.0 and [Ni/H] = 0.0. The tests were made for 2 signal to noise ratios (SNR): 30 and 100 (200 
reference spectra for each SNR) . The values of the dispersions at 1 — cr of the residuals (difference between the 
estimated and the real value) for each parameter are as follows: 





SNR = 30 


SNR = 100 


Teff (K) 


31 


9 


log(g) 


0.14 


0.05 


[Fe/H] (dex) 


0.04 


0.0013 


[Ca/H] (dex) 


0.03 


0.009 


[Ni/H] (dex) 


0.05 


0.017 



The mean results of these Monte-Carlo runs for each parameter show no bias. The dispersions are accept- 
able. Actually, for the Teff determination using the Ha method, the dispersion is at SNR of 100 (respectivel y 
30) about 10 times (respectively 3 times) smaller than the systematic error (estimated by Cavrel et al.l ( 201ll )) 
linked to the physics behind the models used. We note that, of course, the Ha line is not always available for 
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use: one reason is that it can simply not be in the spectral domain used, another is that this Teff determination 
method cannot be used for all stars (cool stars for example). The excitation equilibrium method is then used. 

3 Future work 

The future work to be done on the software is : 

- On the fly reference grid calculation: dynamic call of the SYNTHE software for calculating the reference grid 
directly from SPADES (completed at the time of writing the proceedings) 

- Fix a method for determining micro-turbulcncc 

- Determine the external errors (as opposed to the internal errors determined by Monte-Carlo). One of the 
methods will be the test on known stars (e.g Sun) 

- In its first vcirsion, the software will be fine-tuned to analyze medium to high resolution GIRAFFE spectra of 
Thick Disk stars . 
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